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FOREWORD 


The  Operations  Research  Center  at  the  Massachusetts  Institute  of 
Technology  is  an  interdepartmental  activity  devoted  to  graduate  educa¬ 
tion  and  research  in  the  field  of  operations  research.  The  work  of  the 
Center  is  supported,  in  part,  by  government  contracts  and  grants.  The 
work  reported  herein  was  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-75-C-0556. 

Richard  C.  Larson 
Jeremy  F.  Shapiro 

Co-Directors 

ABSTRACT 

Disaggregation  and  resource  allocation  problems  can  be  often 
formulated  as  convex  knapsack  problems  with  bounded  variables.  In 
this  paper,  we  provide  a  recursive  procedure  to  solve  such  problems. 

The  method  differs  from  classical  optimization  algorithms  of  convex 
programming  in  that  it  determines  at  each  iteration  the  optimal  value 
of  at  least  one  variable.  Applications  and  computational  results  are 
presented . 


1.  INTRODUCTION 


In  this  paper  we  present  a  recursive  method  to  solve  the  following 
continuous  knapsack  problem  with  bounded  variables: 

(N) :  z  “  min  Z  G  (x . ) 

JeJ®  J  J 


Z  x 
JeJ®  J 

-  P* 

(1) 

*bJ  1  XJ 

-ub1 

JeJ* 

(2) 

Vxj 

JeJ® 

where  (• )  for  JeJ®  are  differentiable  convex  functions  on  the  open  convex 

sets  Xj  C  R,  jeJc  respectively. 

Problem  (*5)  arises  in  many  different  settings.  We  briefly  illustrate 
some  of  the  major  applications. 

In  production  planning  and  scheduling,  problem  (II)  arises  in  the 
following  way.  The  allocation  of  production  resources,  as  for  example 
labor  hours  or  machine  capacity,  is  made  with  aggregate  data.  The  results 
are  aggregate  plans  that  must  be  implemented  at  the  detailed  level  and  therefore 
need  to  be  disaggregated.  It  is  at  this  point  that  problem  (N)  plays  an 
important  role.  Two  of  the  objective  functions  encountered  in  these 
situations  are 


S.D.  ^ 
a)  G.  (x, )  -  -tr^  +  4 


*  Vj 


JeJ' 


where  for  each  iten  1,  Sj ,  Dj ,  h ^ ,  and  Xj  correspond  to  the  setup  cost, 
the  demand,  the  holding  cost  per  period  and  the  number  of  hours  of  labor  to 
be  assigned  to  the  item. 


b>  w  -  ( 
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P*  -  *  x  -  AI 

keJ®  _ _ _  _  xj  1 


Z  D. 
keJ®  k 


r*) 


JeJ* 


The  quantity  P°  represents  the  total  number  of  labor  hours  assigned  by  the 
aggregate  plan  to  the  family  composed  by  the  set  of  items  {k:  keJ®}.  For  any 
generic  item  k£J°,  AI^,  SS^,  D^,  and  x^  are  respectively,  the  available 
inventory,  the  safety  stock,  the  demand,  and  the  number  of  labor  hours  to 
be  allocated  to  that  item.  In  this  case,  problem  (N)  attempts  to  equalise  the 
run  out  time  of  each  item  j,  i.e.,  the  length  of  time  that  the  number  of 
units  produced  during  will  last,  with  the  run  out  time  of  the  family 
(k:  k£J°} . 

Constraint  (1)  typically  assures  that  the  disaggregation  process  is 
consistent  with  the  aggregate  plan  while,  the  lower  and  upper  bounds  (2) 
guarantee  respectively  that  the  demand  for  each  item  will  be  met  and  that 
the  overstock  limit  will  not  be  exceeded.  Other  applications  to  production 
planning  can  be  found  in  [2]  and  [3].  A  related  application  is  the  extensive 
use  of  problem  (N)  in  inventory  control  to  dtermine  lot  sizes  and  buffer 
stocks  under  aggregate  constraints  (see  [6]  and  [8]) . 

Very  often  the  allocation  of  financial  resources  gives  rise  to  knapsack 
problems  with  bounded  variables.  A  typical  objective  function  in  this 
instance  is 


G  (x, )  -  -s  (x.+c,)/(x,+m  )  with  m  >  c. 

JJ  jjjjj  jj 


jeJ‘ 


-Gj(Xj)  represents  the  return  of  investing  x^  units  in  activity  J.  Constraint 
(1)  indicates  that  the  total  amount  available  for  investment  must  eaual  P° 
and  the  lower  and  upper  bounds  are  imposed  by  the  market.  Of  course,  the 
function  Gj(Xj)  need  not  be  identical  for  all  j£J°. 

Some  other  contexts  where  problem  (N)  has  been  applied  are  the  theory 


iHtfiSrtrit  -a 
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of  search  [4],  the  allocation  of  promotional  resources  among  competing 
activities  [9],  [11],  and  subgradient  optimization  [5]. 

In  many  instances  where  (N)  plays  an  important  role  the  set  J°  may 
have  several  hundreds  or  thousands  of  elements.  Moreover,  the  problem  has 
to  be  solved  many  times.  In  hierarchical  production  planning,  for  exam? le , 

(see  [2]  and  [3]),  problem  (N)  typically  is  used  in  two  or  more  levels  of 
disaggregation  and  is  solved  at  the  beginning  of  each  period.  Consequently,  an 
efficient  method  to  solve  the  continuous  knapsack  problem  with  bounded 
variables  is  of  central  interest  to  many  applications. 

Formerly,  this  problem  has  been  treated,  for  particular  objective 
functions,  by  convex  programming  arguments  [4]  and  by  dynamic  programming 
[11].  Luss  and  Gupta  [9]  have  presented  an  iterative  method  mainly  for 
strictly  convex  decreasing  functions  and  a  one  pass  algorithm  for  a  set  of 
particular  functions  with  the  variables  hounded  from  below.  Luss  and 
Gupta's  method  consists  in  relaxing  the  upper  bounding  constraints.  The 
algorithm  proposed  in  this  paper  applies  for  a  more  general  set  of  functions 
whenever  the  simpler  problem  of  the  form  min  {£  ’  st  ^xj  “  p»  x  j eXj ) » 

obtained  by  relaxing  both  bounds  (2)  and  solved  at  each  iteration,  has  J*  solu¬ 
tion.  This  condition  is  met  in  most  practical  applications.  We  also  provide 
results  for  the  case  where  the  relaxed  problem  has  no  solution.  At  the  end 
of  each  iteration  of  our  algorithm  (except  possibly  the  last) ,we  show  that 
either  the  subset  of  variables  (xj  :  x^  >  ubj)  have  optimal  value  ubj  in 
(N)  or  the  subset  {xj  :  Xj  £  Slbj}  have  optimal  value  ibj  in  (N) .  This 
characteristic  of  the  algorithm  differs  from  classical  convex  prograrr'inp 
methods.  More  recently  Armstrong,  Cook,  and  Palacios-Gomez  [1]  developed  a 
branch  and  bound  algorithm  to  solve  a  problem  related  to  (N) .  They  require 
the  variables  x^  to  be  either  equal  to  zero  or  to  be  in  the  interval  [l,u^]jeJ* 
and  the  functions  (  )  jeJ®  are  assumed  to  be  concave  increasing  on  [l,Uj] 


and  zero  elsewhere. 


2.  THE  ALGORITHM 


In  section  1,  we  assumed  that  the  functions  (x^ )  are  convex  and 
differentiable.  Without  loss  of  generality  we  can  add  the  conditions: 

a)  ubj  >  fcbj  jeJ®  since,  if  ub^  “  fb^  for  some  keJ®  the  value  of  is 
determined.  Hence,  k  can  be  deleted  from  J#  and  P®  replaced  by  P'-x.  . 

b)  E  Lb,  <  P®  <  I  ub , .  Otherwise  the  problem  is  either  infeasible  or 

jeJ®  J  jeJ®  J 

the  solution  is  trivially  determined. 

(N)  is  a  convex  problem  with  linear  constraints.  Hence,  the  correspond¬ 
ing  Kuhn-Tucker  conditions  (3)  -  (9),  below,  are  necessary  and  sufficient  [10] 


for  optimality  of  jeJ® 


Let  DG^(Xj)  denote  the  derivative  of  G^(*)  at  the  point  x^  and  let 
x*  for  jeJ®  be  an  optimal  solution  to  (N) .  The  Kuhn-Tucker  conditions  are: 


DVx*)  +  i  +  uj  -Tj  -  0 

VV“V  ■ 0 

wv  ■ 0 

E  x*  -  Pe 
jeJ®  J 

“i  i  -J  i 

AeR,  Uj  >  0,  >_  0 

XJ  eXj 


jeJ® 

(?) 

jeJ® 

(4) 

jeJ® 

(5) 

(‘> 

jeJ® 

(7) 

,l£j® 

(") 

jeJ®  • 

(9) 

A,  Uj ,  and  Tj  are  the  Kuhn-Tucker  multipliers  associated  with  the  constraints 
(1),  ubj-Xj  >.  0  and  Xj-fbj  >  0  jeJ®  respectively.  When  ub^  *  +»  [£b^  “  -»] 
for  some  jeJ®,  the  corresponding  condition  (4)  [(5)]  and  multiplier  u^  U T ^ ] 


are  deleted 


We  first  state  and  prove  the  optimality  of  the  algorithm  under  the 
assumption  that  all  relaxations  N(8)  of  (N)  encountered  in  step  1  have  an 
optimal  solution.  At  the  end  of  this  section  we  relax  this  assumption. 


The  Algorithm 

Initialization:  J1  =  J° ,  P1  »  P8. 

Iteration  8  (8*1,2, 3...) 

Step  1:  Solve  N(8)  :  min  {  E  „  G  (x.)  :  E a  x.  ■  P^,  x.eX.  jeJ^} 

jeJ°  3  3  jeJ°  3  3  3 

and  let  the  solution  be  x®  jeJ^.  If  f.bj  £  x^  £  ub^  jeJ^  define 

Xj  =»  xp  jeJp  and  stop, the  solution  x^  jeJ8  generated  by  the 

algorithm  is  optimal.  Otherwise  go  to  step  2. 

Step  2:  Compute 

,8 

V-h^-uDj/  wuece  -  ijt. 

'+ 


*  E  (x^-ub. )  where  *  {jeJ^ 
JeJ®  3  3  + 


-5  i  -v 


and 

nB 


E  o  (f,h,-x®)  where  =  {jeJ^ 
jeJf  J  J 


*i  1 V 


Step  3:  If  £  vB  define  x^  *  ub^  jeJp  and  let 

,8+1 


tB 


Pe+1  -  P6- 


E»  ub, . 

Jejf  3 


If  A^  <  define  Xj  *  Hbj  jeJ®  and  let 


,8+1 


jM.  P6+l 


-  I*  lb.. 
jeJe  3 


0^1  ft 

If  *  4>  stop.  The  solution  Xj  jeJ8  generated  by  the  algorithm 

is  optimal.  Otherwise  let  8  ■  8+1  and  go  to  step  1. 


O 

Since  at  each  iteration  the  set  J  is  reduced  by  at  least  one  element,  the 
algorithm  is  finite.  The  algorithm  relies  on  the  fact  that  it  is  generally 
much  easier  to  solve  N(8)  than  (N).  Problem  N(8)  can  be  solved  by  using 
its  Kuhn-Tucker  conditions.  In  fact,  in  several  practical  applications  it 
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is  even  possible  to  obtain  a  solution  in  closed  form. 

To  prove  that  x*  jeJ8,  generated  by  the  algorithm,  solve  (N) ,  we  construct 
a  corresponding  Kuhn-Tucker  vector  through  the  following  results.  Let  X  be 

the  Kuhn-Tucker  multiplier  associated  to  the  knapsack  constraint  ir.  *’(?). 

8  8 

Lemma  1;  If  at  iteration  6  iP  >  Vp  then 

8  8  8  B 

a)  foranyseJ?  we  have  -DG  (ub  )  >  -DG .  (ub. )  for  all  ieJ-J  and  -DG  (ub  )  >  A 

W  X6+1  <  X* 

Proof; 

a)  A8  “  -DGj (Xj )  jeJ8.  Let  seJ8  and  ieJ8-J8.  Since  the  functions  G^(*) 
are  convex,  for  any  pair  Xj ,  xj  in  we  have  [10]; 

[DGj (x*)  -  DG  <x*)](x*  -  xj)  >  0  (10) 

It  follows  that 


-DGg(ubg)  >  -DGg(x8)  -  A6  ->  -DG1(x^)  >  -DG^ub^ 


b) 


jeJ8-J8 


P*  -  l  x1 

i&l 


8  1 P8  - 


£  ub. 


,6+1 


.6+1 


jeJ 


e+i 


8+1  8+1  8 

For  at  least  one  j  eJp  we  have  x,  >  xp  .  Thus 

O  J  —  1 

Jo  Jo 

A 8+1  =  -DG  (x8+1)  <  -DG  (x8  )  -  A8 

Jo  Jo  ^o  3o 

where  the  Inequality  follows  from  (10). 

Similarly , 

8  8 

Lemma  2;  If  at  iteration  B  Ap  <  Vp  then 

a)  for  any  seJ8  we  have  -DG  (lb  )  <  -DG  (£b.)  for  all  ieJp-Jp  and  -DG  (£b  )  <.  A' 

b)  A8+1  >  A8 

The  proofs  of  this  lemma  and  of  theorem  4  are  omitted  because  they  are  similar 
to  those  of  lemma  2  and  theorem  3  respectively. 


6_t6 
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Theorem  3:  Assume  that  A8  V®,  A*  <  V*  i-0+1,0+2, . . . ,y-l.  Then 

a)  JY  3  (J8-J8-J8) 

b)  X8  >  XY 


Proof: 


»e+i  , 


8+1 


JeJ 


8+1 


JeJ8-J8 


x8  -  E  ub,  +  E  x? 

3  g  J 


J 


jeJ 


6 


jeJ 


E  x8+A8 

j£j6+1 

(11) 


r.8+2 


8+2 


jeJ 


8+2 


E  x1 


6+1 


Jib , 


jeJ 


6+1 


jeJ 


8+1 


E  x8  +  E  x8  +  A8  - 

j  eJ6+1- J8+1  J  jeJ8+1  j  jeJ8+1 


E  tb  . 


.6+2  t8+1  t8+1 

ias,  since  J  *  J  -J 


,3+2 


E  x8  +  A8  -  E  (Jlb.-x8). 


jeJ 


6+2 


jeJ 


6+1 


j  3' 


(1?) 


Similarly  we  obtain 


E  x8  +  A8  -  V  E  (£b  -x8). 
jeJY  j  s‘8+1jEJs  ^  3 


(13) 


From  (11)  it  follows  that  for  at  least  one  j  eJ 

o 


8+1 


8+1  .  6 
xj  1  XJ 

Jo  Jo 


(1*) 


But,  from  the  Kuhn-Tucker  conditions  for  N(8)  and  N(6+l): 


-DG^x8)  -  X8  and  -DG1(x8+1)  -  X84-1  ieJ8+1 


(15) 


Combining  (10),  (14),  (15),  and  lemma  lb)  It  follows  that 

t6+1 


x8*1  1  x8  for  all  jeJ8 


(16) 


Note  that  if  A8  “  0,  then  x8+1  “  xj  solve  N(8+l)  and  if  A®  >  0,  at 

least  one  inequality  in  (16)  will  be  a  strict  inequality.  Moreover,  if  the 

B+l 

functions  G,  are  not  all  strictly  convex,  it  is  possible  that  -DG,  (x,  )  ■ 

J  ■*<>  ^o 

X8+1  -  X®.  In  this  case  N(8+l)  may  have  more  than  one  optimal  solution. 

However,  at  least  one  will  satisfy  (16).  Thus 

J6+1CJ6  and  A6-  Z  (Ab.-x8)>0.  (17). 

'  *  '  Jejf1  3  3 

OiO  D  L-I  ft+1 

Similarly,  (12),  (17),  and  the  fact  that  JP  •  J  -J_  imply  that 

x8+2  1  XJ  all  jeJ8+2,  J®+2  C  J?  and 

A8  -  Z  (Ab.-x8)  -  Z  (Ah  -x8)  >  0. 
jeJ8+1  j  j  j£J8+2  j  j 


Thus, 


Z  xf3  >  Z  x8 
Jet6*3  3  Jd6*3  3 


Continuing  with  the  same  reasoning,  we  obtain 


JiCJ0  i“P+l,...,Y-l;  A8  -  Z  Z  (Ab  -xj)  >  0 
- o  =  «4-1  o  J  J 


3=8+1  jeJ8 


and  from  (13) 

-  2  *}>  2  *5-  (i?) 

jeJY  jeJY 

8+1  8  8  1+1  1  1 

These  conclusions  together  with  J  =  J  -  and  J  “  J  i“8+l, . . . ,Y~1 

prove  part  a).  From  (18)  it  follows  that  xY  ^  x8  jeJY.  Thus,  from  (lft) 
and  (15)  with  Y  instead  of  8+1 


XY  -  -DG^xJ)  <  -DGj(x8)  -  X8  . 
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Theoren  4:  Assume  that  A®  <  V8,  A1  >  V1  l-P+1,0+2, . . . ,y-l.  Then 

a)  JY  2 (J6-J8-J8) 

b)  X8  <  XY 

* 

Theorem  5:  The  set  x^  jeJ  generated  by  the  algorithm  is  optimal  in  (N) . 

Proof:  By  lemmas  1  and  2,  theorems  3b)  and  4b),  the  set  x*  jeJ°  generated 
by  the  algorithm  has  the  following  property: 

-K  (ub  >  »  ...  »  -DGk  (ubk  )  >  -M  <«* >  '  ••• 

11  P  P  11 

-  -DGv  (x*  )  >  -DGi  (Ab±  )>...>  -DG±  (2bt  )  (19) 

g  g  ll  ll'  s  s 

where,  the  indices  k^,  kj .  kp;  i^,  is,  and  v^,  v2,.*.»  v^  correspond 

to  variables  for  which  the  optimal  value  was  set  at  the  upper  bound  and  lower 
bound  in  step  3  of  the  algorithm  and  variables  whose  optimal  value  was 
obtained  in  step  1  of  the  last  iteration  of  the  algorithm  respectively. 

To  see  that  conditions  (3)  -  (9)  are  satisfied  take. 


X  -  -DG  (x  ), 

V  V1 


T  -  X  +  DG  (Jib  )  >  0,  u,  =0  j-1,2, . . .  ,8 

3  3  *j  j 


il  -  -DG  (ub  )  -  X  >  0,  T  -  0  j-1,2 . .  and 

3  J  3  3 


*  uv  “  0  j“l>2 , . . . ,g. 

3  3 


Thus,  since  x^  jeJ°  also  satisfies  (6),  (7),  and  (9)  it  follows  that 

x.  -  ub  j-1,2 . .  x  -x*  j-1,2,. ..,g;  and 

*J  *J  Vj  Vj 


x.  ■  &b  j-l,2,...,s  solve  (N). 

J  J 
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Slnce  the  algorithm  depends  strongly  on  the  existence  of  a  solution 
to  N(B)  we  prove  the  following  theorem.  Let  x^ (A)  denote  a  point  In  P. 
where  DGj[x^(A)]  *  A. 

g 

Theorem  6:  If  (•)  jeJ  are  strictly  convex  and  differentiable  on  X^ 
then,  a  necessary  and  sufficient  condition  for  N(B)  to  have  a  solution  is 
that  there  exist  A^  and  A^  such  that 

l  x,(A,)<Pe<  E  x, (A,)  with 
jeJ°  1  1  jeJ*  J 

Xj(A1)eXj  and  Xj(A2)cXj  jeJB 


g 

Proof :  Sufficient  Condition;  The  functions  (A)  jeJp  are  continuous. 

g 

Thus,  Z  _x, (A)  is  a  continuous  function  of  A.  Moreover,  X.  jeJP  are  open 
Tp  J  J 


jeJ 


intervals  in  R.  Therefore,  there  are  a  AqeR  and  points  Xj(Ao)eX^  je.T  such 

that  Z„x.(A)»PB. 
jeJB  i 

Necessary  Condition:  Follows  from  the  Kuhn- Tucker  conditions 
for  N(B) .  (H 

It  is  worth  pointing  out  that  a  necessary  condition  for  N(B)  to  have  a 
solution  is  that 


max  lim 

j«e  V*j 

*jeXj 


"VV 


<  min  lim 

)€J6  *j'‘j 

Vxj 


"VV 


(2*) 


where 

g 

Sj  »  inf{xeR:  x^eX^}  JeJp  and 

g 

bj  -  sup{xcR:  XjEXj}  jeJP 

Condition  (20)  can  be  useful  as  a  tool  to  conclude  that  N(B)  has  no  solution. 
When  there  are  strictly  increasing  and  decreasing  functions  over  X^  among 

g 

the  functions  G^(*)  jeJ  ,  N(B)  has  no  solution.  In  this  case  the  left  band 
side  of  inequality  (20)  is  negative.  The  next  theorem  shows  how  to  cope  with 


n'-rmhm  waaSMeiSiefl. 
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this  difficulty. 


Let  ■  {jeJ°s  Gj(*)  is  strictly  increasing},  0  0  and  let 

J2  =  tjeJ®:  Gj(*)  is  strictly  decreasing},  J2  0  0.  Note  that  we  do  not 
require  J°-J^-J2  to  be  an  empty  set.  Assume  that  (N)  has  an  optima’  solution 
Xj  jeJ®,  ubj  <  +00  and  ibj  >  -°°  jeJ®. 

Theorem  7:  The  optimal  solution  of  (N)  is  such  that  either 

x =  ub^  ieJ2  and/or  x^  ■  ib^  ieJ^. 

Proof :  The  optimal  solution  satisfies  (3)  -  (9)  and  in  particular 

-DG^  (ubj )  -  X  +  Uj  j  eJ  (ub)  =  { j  eJ® :  x^  ■  ub  ^  } 

-DGjCfcbj)  -  X  -  Tj  jeJ(ib)  -  {j eJ° :  x*  -  *b  }  and 

-DGj (x*)  -  X  jeJ®  -  J(£b)  -  J(ub). 

Note  that  the  assumption  P®  >  £  £b.  implies  J®~J(£b)  0  0.  If 

*  j 

min{-DG^ (Xj)  jeJ°-J(db)}  0,  since  -DG^tx^)  <  0  ieJ^,  it  follows  that 
CJ(tb).  Otherwise,  since  -DG^tx^)  >  0  ieJ2  and  -DGj(ub^)  >_  -PG^Cx^} 
-DG^CAb^)  for  any  jeJ(ub)  ieJ*-J(ub)-J(ib)  and  keJ(fcb),  we  have  that 
J2  Cj(ub).  M 

A  direct  consequence  of  theorem  7  is  that  the  solution  to  (N)  car.  be 
obtained  by  solving  the  following  two  problems 

(N  ):  z  -  mint  £  G,  (x.):  £  x  =  P°-  £  ub  S-b  <x  <ub,  JeJ®-J,} 

1  1  jeJ®-J2  3  3  jeJ®-J2  3  jeJ2  J  J  J  J  2 

(N,):  z„  ■  mint  £  G, (x,):  £  x.  *  P®-  £  £b,  Jtb  <r 4<v>*  jeJ°-J-} 

2  3  3  jeJ’-Jj^  3  jeJx  3331  1 

and  taking  z  ■  min  (  £  G  (ub,)+z  ,  £  G  (£b. )  +  z«).  By  analyzing  tbe 

jeJ2  3  3  jeJj  3  3 

Kuhn-Tucker  conditions  (3)  -  (9)  it  becomes  apparent  that  problem  (N)  has 
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no  solution  if  there  are  simultaneously  a  strictly  Increasing  G.  (•)  with 

31 

£b,  -  and  a  strictly  decreasing  G.  (•)  with  ub.  -  +°°  and  these  two  limits 

31  32  32 
are  attainable  on  X.  and  X  respectively.  Our  algorithm  does  not  apply  in 
J1  32 

these  instances.  However,  it  is  not  difficult  to  show  that  Theorem  7  holds 
if  only  decreasing  functions  (increasing  functions)  have  unbounded  upper 
bounds  (lower  bounds). 

Concluding  this  section  we  consider  a  version  of  problem  (N)  when  the  knapsack 
constraint  is  an  Inequality.  Let  (N<)  and  (N>)  be  the  versions  of  (?!)  with 
constraint  1  being  an  inequality  of  the  type  £  and  ^  respectively.  Tor  each 
index  jeJ4  let  h^  be  the  value  of  Xj  for  which  DGj(xj)  -  0  over  X^.  If  such 
point  does  not  exist  we  adopt  hj  -  -*°  (+")  in  Theorem  8  (Theorem  9)  below. 

K 

Let  Xj  jeJ*  solve  (N).  The  indices  k^ ,  v^ ,  and  i^  correspond  to  those  in 
expression  (19). 

Theorem  8: 

a)  If  A  «  -DG  (x*  )  >  0,  x,  *  x*  JeJ4  solve  (N<) 

VJ  3  33 

b7  If  »  -DG  (x*  )  <  0,  x,  defined  as: 

VJ  V3  3 

x.  "  j-l,2|...|S; 

3  3 

*  "  max  (Lb  ,h  )  3“l»2,...,g; 

j  3  v3 

x,  •  max  (£b.  ,h.  )  for  all  k,  such  that  DG.  (ub.  )  >  0;  and 
j  j  i  3  3  3 

x.  ■  ub.  for  all  k,  such  that  DG.  (ub.  )  <  0  solve  (N£) . 

j  3  3  K3  3 

Theorem  9: 

a)  If  A  -  -DG  (x*  )  <  0,  x  -  x?  jeJ9  solve  (N^) 

v3  v3  3  3 

b)  If  -  -DG  (x*  )  >  0,  x  defined  as: 

V3  3  3 

x  -  ub  3-1, 2,..., p; 

j  3 


X  “  mln(ub  ,h  )  j-l,2,...,g; 

Vj  Vj  Vj 

x.  “  min(ub  ,h  )  for  all  1.  such  that  DG  (f-b  )  £  0;  and 

i  3  i 

x.  “  lb  for  all  1  such  that  DG  (lb  )  >  0  solve  (N^). 

*J  3  i  i 

Theorems  8  and  9  show  that  problems  (N<)  and  (N>)  can  be  solved  by  first 
applying  our  algorithm  to  solve  (N)  and  next  making  appropriate  chanpes  to 
x*  jeJ.  The  proofs  of  both  theorems  have  been  omitted  because  of  their  simpli¬ 
city .  The  extensions  of  the  last  two  results  for  the  case  where  the  solution 


to  problem  (N)  has  all  variables  at  the  upper  or  lower  bound  Is  strnlpht^orward. 
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3.  COMPUTATIONAL  RESULTS 

Tables  1  through  7  present  the  results  of  solving  84  problems  of  type 
(N) .  All  problems  were  solved  by  the  algorithm  described  In  this  paper  and 
whenever  applicable  by  Luss  and  Gupta's  method  [9].  For  Identification 
purposes,  our  algorithm  Is  denoted  by  BH,  while  the  algorithm  in  [9]  is 
denoted  by  LG.  In  each  problem,  the  objective  function  was  composed  by 
functions  G^(x^)  of  the  same  functional  form.  They  are  indicated  in  the  first  row 
of  each  table.  Following  the  time  in  seconds  to  solve  each  problem,  by  both 
methods,  is  the  number  of  iterations  required,  n  represents  the  number 
of  variables  in  a  problem.  For  a  fixed  n  we  solved  three  problems  for 
each  type  of  objective  function.  In  Luss  and  Gupta's  method  [9]  the  order¬ 
ing  of  the  derivatives  evaluated  at  the  lower  bound  of  each  variable  was 
executed  by  the  "Quicksort  Method"  [7].  In  our  algorithm,  problems  T.T(B) 
were  solved  using  the  corresponding  Kuhn-Tucker  conditions.  Luss  and  Gupta’s 
algorithm  does  not  apply  to  the  problems  of  Tables  6  and  7,  because  either 
the  objective  functions  are  strictly  convex  increasing  (Table  6)  or 
we  have  not  imposed,  as  their  algorithm  requires,  any  condition  among  the 
values  of  the  bounds,  P°  and  the  point  where  each  of  the  G^(x^)  attains  its 
minimum  (Table  7) .  The  computer  used  is  a  Borroughs  B6700.  The  programs  were 
written  in  Algol.  Applications  of  problems  presented  in  Tables  5  and  7  to 
hierarchical  production  planning  can  be  found  in  [2]  and  [3].  The  para¬ 
meters  (s^.m^.Abj.ubj,  etc.)  corresponding  to  problems  of  Tables  1  through 
5  (6  and  7)  were  randomly  generated  in  intervals  where  the  functions  ^(x^) 
are  strictly  convex  decreasing  (strictly  convex). 

For  the  problems  presented  in  Tables  1  through  5  we  have  noticed  that  the 
time  required  by  Luss  and  Gupta's  method  just  to  compute  the  derivatives  at 
the  lower  bound  and  order  them  la  comparable  to  the  total  time  of  our  algorithm. 


As  a  final  note,  we  would  like  to  point  out  that,  although  we  have 
tried  to  program  both  algorithms  as  efficiently  as  we  could,  the  computational 
results  should  be  looked  upon  with  caution.  The  intent  of  including  the 
results  of  these  experiments  is  to  give  to  the  reader  an  idea  of  the  tine 
required  to  solve  problems  of  the  type  (N)  with  different  objective  functions 
and  sizes.  The  major  attraction  of  our  algorithm  is  its  applicability  to  a 
wider  class  of  problems  than  the  one  in  {9]. 
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